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Abstract 

Using two-color lattice QCD with Wilson fermions, we report a study of the finite 
baryon number density system with two-flavors. First we investigate the Polyakov 
line and thermodynamical quantities in the (k, fj,) plane, where K and /i are the 
hopping parameter and chemical potential in the fermion action, respectively. Then 
we calculate propagators of meson (qTq) and baryon (qTq) states. We find that the 
vector meson propagators are strongly modified in large jx regions, indicating the 
reduction of the mass. This anomalous behavior of the vector meson is observed for 
the first time in lattice QCD. 

PACS: 12.38Gc, 11.15Ha, 12.38Mh, 24.85+p 
Keywors: Lattice gauge theory, Finite density QCD 



1 Introduction 

QCD has confinement and deconfinement phases when the temperature, T, is varied. 
When an additional parameter, the density or the chemical potential, is added, QCD 
may have a much richer structure |l|, ||] . Experimentally, a wide region of the (T, //) 
plane has been investigated by AGS, SPS and RHIC, and a higher density realm 
might be covered by GSI and JHF in the future. There have been experiments 
which suggest that vector meson masses are modified at finite density Q Q . 

Lattice QCD is expected to provide nonperturbative information of the hadronic 
world at finite density as the first principle calculation ||. However, numerical study 
of lattice QCD with chemical potential is extremely difficult, because at finite fj,, 
the fermion matrix does not satisfy the usual condition, 

D ] = 75^75, (1) 



1 



and the fermion determinant det-D becomes complex, which appears in the Eu- 
clidean path integral measure, 



Z = J P[/PWe-^-* D * = J VUdetDe~P Sa . (2) 

In order to circumvent the above difficulty, two-color QCD has been investigated 
@, 0, |, |, 10, 11, [12| . Although recent progress of lattice calculations with the finite 
chemical potential has been prominent ||, 0, ||, |l|, , it is still very difficult to 
study regions around critical (i at low temperature by lattice QCD simulations. We 
report here, for the first time, a lattice study of hadron propagators together with 
thermodynamical quantities for finite density SU(2) QCD with Wilson fermions. 
Two-flavor case is studied. 

The two-color SU(2) theory qualitatively has most of the important features 
of QCD, such as deconfinement transition at finite temperature. The essential 
difference between two- and three-color systems is that in the SU(2) case baryons 
are made of two quarks, i.e., they are bosons. Therefore we must be careful in the 
study of two-color QCD in order to gain some insight into a subset of the phenomena 
expected in QCD. It is indispensable to compare results for mesons and diquarks 
which have a special relation in the SU(2) case. Recent theoretical progress has 



greatly improved our understanding of the two-color QCD system [fLq, 19, 20, \l 



2 Actions and algorithm 

We introduce the chemical potential in the conventional manner, 
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D(x,x') = 6 XjX , - kJ2 {(1 - li)Ui(x)5 x , >x+} + (1 + 7i)U}(x')6 x , x _ l \ 

i=l 

-« - l±)Ui{x)8 xljX+i + e-"(l + T^V^-i} ■ (3) 

For the gauge action, we employ the plaquette and Iwasaki improved actions. Here, 
we report the improved action case only, and for analyses with the plaquette gauge 
action, we refer to Ref. p2[| . Little is known about dynamical fermion simulations 
in which the chemical potential is introduced. We therefore employ an algorithm 
where the ratio of the determinant, 

is evaluated explicitly in each Metropolis update process, C7 — > C7 + AU, where 



AD = D(U + AU) - D{U). For details of the algorithm, see Refs.g |2J. This 



algorithm has a long Markov step and is very reliable [25|. Numerical costs are, 
however, huge and we are restricted to small lattices. In the following studies, 
therefore, we check that the results obtained are not sensitive to the boundary 
conditions. 
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3 Study of the (ac, ji) Parameter Space 



Since there are few color SU(2) lattice studies using Wilson fermions with finite fi in 
the literature except Ref.||, we first investigate the relevant parameter space. We 
measure the Polyakov line, (L) on a 4 lattice by changing (3 for \i = and k = 0.150 
and choose the region where {L) is small, i.e., the system is in the confinement phase 
at zero baryon number density. We set (5 = 0.7 on the basis of this analysis. 

At this value of (3, we measure (L), its susceptibility, g, the gluon energy 

density, (E g ) =< a(i/T) ^ G ^ an< ^ *he numDer density, (rt) = y-J^\ogZ, as a 
function of [i and k. Here Sq is the gauge action and V s = N x N y N z is the spatial 
volume of the lattice. In Fig.[E], we show (L), (n) and (E g ) on a 4 4 lattice as a 
function of \x and k. They increase as fi becomes large and show the deconfinement 
behavior. We observe that the simulation always breaks down when we increase /i 
further. 




Figure 1: The Polyakov line (L), the number density (n) and the gluon energy (E g ) as 
a function of k and /x. Lattice size is 4 4 . 

In Fig.|2[ we show the Polyakov line susceptibilities k = 0.150 and 0.175 under 
the antiperiodic spatial boundary condition and for k = 0.160 under the periodic 
boundary condition. All exhibit a peak when [i increases, which indicates a decon- 
finement transition. All quantities support the picture that at large fi the system 
undergoes the transition from the confinement to the deconfinement phase. In ad- 
dition to the increase of (L) and (E g ), the rapid increase of (n) is observed. The 
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Figure 2: Polyakov line susceptibility, d(L)/d/fi — {{L— (L))(n— (n))), as a function of /x 
on 4 3 x 8 lattice, (a) k = 0.150, 0.175 under the anti-periodic spatial boundary condition 
and (b) k = 0.160 under the periodic boundary condition. 



instability at large [i may be an indication of a new phase with the diquark conden- 
sation. [ 

Although the behavior of all quantities supports the existence of the deconfine- 
ment phase at large //, there are some indications that suggest a more complicated 
phase. In many cases, we observe a second peak in the Polyakov line susceptibility 
at large \x. 



4 Hadron propagators 

We calculate correlations of color singlet hadron operators, M(x) = rp^ [ (x)T a j3'ipg(x) 
and B{x) = e ab ^(x)(CT) a pfij; b p(x), where T is the product of Dirac matrices and 
f is a Pauli matrix acting on flavor indices. C is the charge conjugation matrix and 
a and b are color indices. We set f = T2 or T2T so that the wave function is totally 
antisymmetric. 

To our knowledge, no study of the behaviors of hadrons including vector mesons 
at finite baryon density has been performed using lattice QCD. Vector mesons are 
important since they provide information at several stages of heavy ion collision in 
the form of lepton pairs. 

In Fig.||, we show propagators of the pseudoscalar and vector mesons and those 
of the scalar (r = 75) and pseudovector (7^) baryons for \i = and \ia = 0.8, 
where a is the lattice spacing. At fj, = 0, the propagator of pseudoscalar mesons 
is equivalent to that of scalar baryons, and that of the vector meson is equivalent 
to that of the pseudovector baryon. We see that this relation is satisfied in the 
numerical calculation. The most prominent feature here is that the vector meson 
propagator is strongly modified at fia = 0.8. Its slope is more gradual than that for 
the pseudoscalar, i.e., the vector meson becomes lighter. 

The lattice size of 4 3 x 8 is too small to extract information of the mass pole. 

1 Since our fermion action includes only bilinear terms of ip and ip and not those of ip and ip, or ip and 
ip, the diquark condensation cannot emerge. 
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Figure 3: Hadron propagators at fia = (left) and /ia 
mesons and scalar and pseudovector baryons are shown. 



0.8. Pseudoscalar and vector 



Nevertheless, we fit the data to clarify the chemical potential effect qualitatively, 
and in Fig.|| we plot the scalar and vector meson masses together with those of cor- 
responding baryons, as a function of the chemical potential. We find that the vector 
meson mass drops as p reaches the critical region. In order to confirm this unex- 
pected result, we calculate both periodic and antiperiodic boundary conditions and 
several k's (k = 0.150,0.160,0.175); the reduction of vector meson mass is always 
observed. Although our lattice size is too small and the statistics are insufficient for 
extrapolation to the chiral limit, the signal of the anomalous behavior of the vector 
channel is clearly seen. 

At p = 0, because of the QCD inequality |Q, the lightest mass should be in 
the pseudoscalar channel. The inequality holds under two assumptions, i.e., (i) 
there is no disconnected diagram and (ii) Eq.([l]) is valid. The second condition is 
not satisfied for the finite density state. Indeed, there are several conjectures in 
the literature. Brown and Rho first proposed the scaling law m*/m = f*/f-K to 
explain the large low mass lepton pair enhancement observed in CERES0, where 
m* and /* are the mass and the pion decay constant in the medium p6[. Based 
on the QCD sum rule, Hatsuda and Lee predicted a decrease of p meson mass as 



a function of p [27]. Harada et al. showed that vector meson mass vanishes at the 
critical density as a consequence of an effective theory with hidden local symmetry 
1 28]. Yokokawa et al. proposed simultaneous softening of a and p mesons associated 
with the chiral restoration [^] . If the sudden drop of the vector meson mass is not 
a special feature of the color SU(2) model, this may be the first lattice QCD result 
to show the reduction of the vector meson mass in the medium. 

We do not observe any special feature in baryon (diquark) channels. We fit the 
baryon propagators as 



G(t) = de-^- 2 ^ + c 2 e- ( - m+2 ^ Nt 



-t) 



(5) 



where C\ and C2 are not independent because of the boundary condition, G(0) = 
G(Nt). At p > 0, the pseudovector baryon propagator is not equivalent to the 
vector meson, but their masses are similar. 
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Figure 4: Meson mass (left) and diquark mass (right) as a function of \i at k — 0.160 
with the periodic boundary condition. 



5 Concluding remarks 

In this study, we have investigated the Polyakov line and its susceptibility together 
with thermodynamical quantities for Wilson fermions in (k, /x) parameter space. 
Since there has been no lattice study using Wilson fermions and improved gauge 
action for the two-color finite density system, this step is necessary in order to un- 
derstand in which region hadron propagators in the medium change their behavior. 

We observed a sudden reduction of the vector meson mass when the chemical 
potential was increased. If this is due to the mixing of meson and diquark states, the 
mass of diquark partner, ip c j^ip, should increase, but this is not the case. Rather, 
ipcl^ behaves in a similar way to the vector. As an urgent project, we will per- 
form simulations on larger lattices, which will allow us to perform the chiral limit 
extrapolation to estimate the physical scale. 

Our lattice here is small, but results for the case of periodic and antiperiodic 
boundary conditions in the spatial directions show the same qualitative behavior. 

The behavior of the thermodynamic quantities together with that of the Polyakov 
line supports the standard picture, i.e., QCD undergoes a transition from the con- 
finement to the deconfinement phase. We have observed several indications in the 
susceptibility and gluon propagators, which may suggest a more complicated phase 
at finite baryon number density. The number density may play an essential role 
in confirming the situation; it should increase from \i = Mb/N c , where Mb is the 
lightest baryon mass and N c is the number of colors. At the deconfinement tran- 
sition point, fj, c , particle which bears the baryonic charge is changed from hadrons 
to quarks, therefore the mass of particles corresponding to the chemical potential 
changes. This results in the jump of the number density at [i c . In the following 
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paper, we will perform detailed study of (n) and its susceptibility. 
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